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ABSTRACT 

We investigate the practice of assigning high spin temperatures to damped Lyman- 
a absorption systems (DL As) not detected in Hi 21 -cm absorption. In particular, 
Kanekar & Chengalur (2003) have attributed the mix of 21-cm detections and non-detections 
in low redshift (z a bs ^ 2.04) DLAs to a mix of spin temperatures, while the non-detections 
at high redshift were attributed to high spin temperatures. Below z a b s = 0.9, where some 
of the DLA host galaxy morphologies are known, we find that 21-cm absorption is normally 
detected towards large radio sources when the absorber is known to be associated with a 
large intermediate (spiral) galaxy. Furthermore, at these redshifts, only one of the six 21-cm 
non-detections has an optical identification and these DLAs tend to lie along the sight-lines 
to the largest background radio continuum sources. For these and many of the high redshift 
DLAs occulting large radio continua, we therefore expect covering factors of less than the 
assumed/estimated value of unity. This would have the effect of introducing a range of spin 
temperatures considerably narrower than the current range of AT S > 9000 K, while still sup- 
porting the hypothesis that the high redshift DLA sample comprises a larger proportion of 
compact galaxies than the low redshift sample. 

Key words: quasars: absorption lines - cosmology: observations - cosmology: early Uni- 
verse - galaxies: ISM 



1 INTRODUCTION 

Despite their relative paucity, the highest column density absorp- 
tion systems, damped Lyman-a systems (where Nm ^ 2 x 10 20 
cm -2 ), are an important component of the high redshift (z ~ 2 — 6) 
Universe since they acc ount for most of the neutral gas available 
for star formation (e.g. lLanzetta et al]|l99ll) . The Hi 21-cm hy- 
perfine transition can provide an alternative and complementary 
view of DLAs. Presuming that the 21-cm and Lyman-a absorption 
arise in the same cloud complexes iDickev & Lockman 1990) 1 , the 
column density [cm -2 ] of the absorbing gas in a homogeneous 
cloud is related to the velocity integrated op tical depth, where 
r = -In (l - -fg), of the 21-cm line via <Wolfe & Burbidgd 
Il975t) : 

N m = 1-823 x 10 18 T spin / V dv , (1) 

where T s [K] is the spin temperature, a is the depth of the line (or 
r.m.s. noise in the case of a non-detection) and S and / the flux 
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1 Where the 21-cm and Lyman-a: absorption trace the cold and total neutral 
hydrogen, respectively. 



density and covering factor of the background continuum source, 
respectively. 

As summarised bv lWolfd fl98ch . 21-cm absorption line stud- 
ies can face some difficulties compared with optical work: (1) in 
blind surveys the bandwidth of radio receiver/backend combina- 
tions only allows small redshift intervals to be searched for H I ab- 
sorption; (2) terrestrial interference can be severe below 1.4 GHz; 
(3) if the spin temperature of the absorbing gas is high or if the 
background source size exceeds that of the absorption cloud(s), 
then the apparent optical depth can be considerably overestimated. 
Despite these technical difficulties, 21-cm studies of DLAs have 
provided information about the physical conditions of absorbers, 
including the kinematics, temperature and velocity distribution of 
the gas. 

One unresolved question concerns the typical size and struc- 
ture of DLAs, with models ranging from large, rapidly rotating 
proto-disks (e.g. Prochaska & Wolfe 111 9971) to small, merg ing sub- 
galactic systems (e.g. Haehnelt, Steinmetz & Rauch 1998). Direct 
imaging by various groups (e.g. iLe Brun et alJll997t iRao et alJ 
2003; IChen & Lanzettall2003h of the low-z DLA host-galaxies re- 
veals a mix of spirals, dwarf and low surface brightness (LSB) 
galaxies. This 'mixed morphology' picture of DLA host-galaxies 
appears to be reflected in the nearby Universe in a recent 21-cm 
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Figure 1. Top: Spin temperature against absorption redshift as given in 
KC03. The solid symbols represent the 21 -cm detections and the shapes 
represent the type of galaxy with which the DLA is associated: circle- 
unknown type, star-spiral, square-dwarf, triangle-LSB (obtained from 
KC03 and references therein). This key is also applied to Figs. l2*lto 151 
The arrows show the spin temperature lower limits and all of these bar one 
(0454+039 at jz a t, s = 0.8596) have unknown host identifications. Bottom: 
The covering factors applied to derive the spin temperatures. The shapes are 
coded as above with the solid and open symbols representing the covering 
factors applied for the 21 -cm detections and non-detections, respectively. 



z = emission study iRvan- Weber. Webster & Stavelev-Smithl 

2003). 

From a survey of 10 sources, in co njunction with the available 
literature. lkanekar & ChengaluJ l2003l hereafter lKC03l) find both 
detections and non-detections of 21 -cm absorption in DLAs at z < 
2, whereas the high-z results are almost exclusively non-detections 
(Fig. Q top). Additionally, many of the z < 1 DLAs have host- 
galaxy identifications from optical imaging, given in Table □ By 
combining the 21 -cm optical depth with the total Hi column den- 
sity from optical spectroscopy covering the damped Lyman-a line, 
KC03 find low spin temperatures (T B ~ 200 K), which are typi - 
cal of the Milky Way and local spirals l Dickev & Lockman 1990), 
in all cases where the absorber was identified as a spiral (Fig. 0. 
Furthermore, all DLAs at z < 1 with temperatures T s > 1 000 K 
were found to be associated with dwarf or LSBs and KC03 inter- 
preted the non-detections at high-z (> 2) as DLAs with high T s . 
They therefore advocate a picture where the DLA host-galaxy pop- 
ulation is dominated by the warmer dwarfs and LSBs at high-z and 
evolves to include a higher proportion of spirals at low-z, thus sat- 
isfying one expectation of hierarchical galaxy formation scenarios. 
However, in many of these studies, particularly at high redshift, the 
covering factor, /, of the background continuum is often assumed 
to be unity (Fig.Q bottom - see also Table Q). Since / plays an 
equal role as T s in relating the column density to the observed op- 
tical depth profile, it is important to consider the robustness of this 
assumption, which is the aim of this article. 



2 FACTORS AFFECTING THE DETECTION OF 21-CM 
ABSORPTION 

2.1 Previous 21-cm absorption searches 

In TableQ we summarise the previously published 21-cm searches 
to which we add the 3 non-detections detailed in Appendix A. The 
optical depths are quoted for various resolutions between 0.6 and 
17 km s -1 . As discussed in detail by Curran et al. (2002a), the 
choice of resolution affects the upper limit to the optical depth ob- 
tained and so we normalise these to 3 km s _1 , a fairly typical res- 
olution for the DLAs detected in 21-cm absorption. In the table 
we also give the covering factors used to determine the spin tem- 
perature in each of the previous searches. Where no high resolu- 
tion radio images are available for both the background source and 
absorber, it is common practice to assume that the absorber com- 
pletely covers the continuum source (e.g. lLane & Briggsl200ll) . al- 
though in some cases the covering factor is estimated as the ra- 
tio of the compa ct unresolved c omponent's flux to the total radio 
flux (e.g. lBriggs & Wolfdl 19831) . A good example of this method 
is given bv lLaneetaTT T2000). where the constancy of the 21-cm 
absorption profile in the z a b s = 0.09 DLA over the partially 
resolved core of 0738 +313 gives an estimate of / ~ 0.98. As 
lBriggs&Woifel<1983l) mention, however, even though this method 
provides information on the background source structure, it pro- 
vides none on the size of the absorbing region 2 . The only way to 
unambiguously determine the covering factor is by mapping both 
the source and the absorbing gas at high angular resolution, which 
is challenging due to the very high sensitivities required. 



2.2 Column density and background flux biases 

In order to check that the non-detection of 21-cm absorption in 
DLAs is not the result of lower atomic hydrogen column densities, 
we plot the velocity integrated optical depth against column density 
in Fig. [2] We use this rather than the peak optical depth as this 
provides a better representation of the strength of the line. For a 
single cloud in thermodynamic equilibrium, the spin temperature 
of the gas can be estimated from the kinetic temperature according 
to T spin w T kin < 2 2 x FW HM 2 (e.g. lLane & Briggj200ll) and 
this has been used bv lKC03| in order to derive spin temperatures and 
line-widths which are self consistent. However, in the absence of 
detailed knowledge of the covering factor, it is impossible to assign 
a specific spin temperature to the system. We therefore assume a 
FWHM of 20 km s" 1 (the mean value of the detections) for the 
non-detections (Table 0. From the table see that the FWHM can 
range from 4 to 50 km s _1 , each of which would not cause too 
large a deviation from our chosen value of 20 km s _1 on the log 
plots. This is confirmed by the fact that using the peak rather than 
velocity integrated optical depths gives similar qualitative results in 
Figs.|2land|3] 

Although the column densities for the non-detections are gen- 
erally lower, they appear to have been searched to correspondingly 
lower optical depth limits. To check this, in Fig.[3]we show the dis- 
tribution of the velocity integrated optical depth modified accord- 
ing to the column density of each DLA: In the optically thin regime 
(<t <C f-S) equationQreduces to 3 



2 The effect of background source size is clearly demonstrated, for i n- 
stance, in absorption studies of PSR B 1 849+00 1 Stanimirovic et al. 2003). 

3 As seen from TablelTI this is a reasonable assumption providing that / is 
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Table 1. The DLAs and sub-DLAs which have been searched for 21 -cm absorption. In the top panel we list the detections and in the bottom panel the 
non-detections. T pe ^ is the peak optical depth of each individually resolved component, the number of which are listed under n comp (B designates 
a blend) and Av is the full width half maximum (FWHM) of the line [km s ], For the non-detections the 3cr upper limits of r pcak at a velocity 
resolution of 3 km s~ 1 (see main text) are quoted. / and T s are the covering factor used and the inferred spin temperature [K] from the literature, 
respectively. For the former + denotes that / is estimated from high resolution radio images and * denotes that / is assumed. TVjji and jz a t, s are the 
total neutral hydrogen column density [cm -2 ] (see http://www.phys.unsw.edu.au/~sjc/dla i and the redshift of the DLA, respectively, with the optical 
identification (ID) given: D-dwarf, L-LSB, S-spiral, U-unknown. z C m is the redshift of the background QSO and S is the flux density [Jy] at ^ b s , 
obtained from either the reference given or as cited in Curran et al. (2002b). The final column gives the 21-cm absorption reference. 
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iBrigp & Wolfelll983h. ffllWolfe et idT{l985h (8llCarilli et al]il992h. (9)lLanzetta et alJll995ft. qO^ICarilli et alJil996h. (llllde Bruvn et alJll996l) . 
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Figure 2. The velocity integrated optical depth versus total neutral hydro- 
gen column density for the DLAs searched for 21 -cm absorption. In this 
and the following figure, the shapes/bold histogram represent the detections 
and the arrows/hatched histogram the upper limits. 
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Figure 3. The modified velocity integrated optical depth versus the radio 
flux density at i/ obs for the optically thin (t 5; 0.3) 21-cm searches. 

N m = 1.823 xl0 18 ^P J~dv, (2) 

and from TableQwe see that this is applicable (r < 0.3) to all but 
two systems, 0235+164 and 0458-020. With this approximation we 
can introduce a linear correction to the optical depth in order to ac- 
count for the range of column densities exhibited by the DLA sam- 
ple. Here we correct the optical depth by a factor of 2 x 10 20 /-/Vhi, 
where Nm ) 2 x 10 20 cm -2 is nominally defined as the limiting 

not small. For / ^ olTs me assum pti° n i- s n °t valid, although it supports 
the argument that / < 1. 



column density of a DLA. This has the effect of putting the data 
points on an equal footing on the ordinate, i.e. for a given optical 
depth limit we account for the fact that a non-detected system of 
low column density has not been searched for as deeply as one of 
higher column density. From Fig. [3] we see, not surprisingly, that 
detections tend to occur towards higher flux densities, although the 
samples are not mutually exclusive. For the modified peak optical 
depths the overlap is more pronounced than in Fig. [3] and for both 
samples —3 < log 10 r + [20.3 — log 10 Nm] ~ —1.5. Equation|2| 
therefore suggests similar values of //T sp i n (over a 2 order of mag- 
nitude range) for the detections and non-detections. We discuss the 
possible contributions of / and T sp i n in the following section. 

2.3 Spin temperatures and covering factors 

From Section l2~2l it appears that neither the column density nor 
background flux (providing they are sufficiently high) are the main 
factors in determining whether a DLA is detected in 21-cm absorp- 
tion. As previously discussed, the popular consensus is the wide 
range of estimated spin temperatures (possibly 20 to > 9000 K) 
in DLAs. In view that the equally influential covering factor is often 
assumed to be unity, we note the following: 

(i) In radio quasar hosts the detection rate for 21-cm absorption 
is muc h higher in the comp act sources, wher e the coverage is high 
IPihlstrom et alj2003llvermeulen et all2003l) . 

(ii) The presence of H2 in 0528-020 

iLevshakov & Varshalovichl 1 1 9851) suggests that the gas i s 
relatively cold with T kin « 200 K Jgrianand & Petitieanlll998h. 
cf. r BP in > 700 K from a 21-cm non-detection ICarilrTeHtl] 
1996). This of course could be the result of different sight-lines 
being probed as well as the possibility the absorption is not due 
to a single cloud. This is known to be a complex situation with 
Zabs ~ Zcm (Table0, suggesting an infalling system rather than a 
simple intervener, although it is one of the cases where a covering 
factor of unity is assumed. This is a valid assumption, however, if 
the absorbing material has an extent of > 70 pc. 4 

Note that a similar situation is found towards 0738+313, where 
Tkin ~ 300 K, cf. r s p in « 800 K in both absorbers lLane et alJ 
2000: iKanekar et al. 2001 |). As seen from the absorption profiles 
towards the 1.3 GHz VLBA image (figure 2 of Lane et alj|200fj|) . 
although overall / is large, the coverage towards the north appears 
to be considerably less than towards the south-east. This suggests 
that the absorber extends to < 16 pc and « 22 pc in these respective 
directions. A lower coverage in the northerly direction could pos- 
sibly reduce the spin temperature estimates to closer to the kinetic 
te mperature values, although the two-phase gas model suggested 
bv lLaneetaD f2000) offers an alternate explanation. In this model, 
the deep, narrow absorption lines give the temperature of the cold 
gas, and the much wider weaker line giving the warm gas tempera- 
ture. These models are discussed further at the end of this section. 

(iii) For the DLAs in which 21-cm absorption has been detected 
there appears to be little correlation between spin temperature and 
line-width. However, these are complex, multi-component systems 
characterised by a single velocity integrated optical depth and con- 
sidering the possibility of both narrow/cold and wide/warm compo- 
nents in the profile (above), the lack of correlation is not surprising. 
Such simplification does however demonstrate the possible pitfalls 
in estimating line-widths for the non-detections (Section f2.2t . 

From the radio continuum size of 0.01" (see TablefSJ at z = 2.8 (H = 
75 km s" 1 Mpc" 1 , fWtef = 0.27). 
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Figure 4. The angular sizes of the background radio continuum sources. Left: The dashed lines represent the angular extent of different size absorption clouds 
(D c ) as a function of redshift (Hq = 75 km s — 1 Mpc~\ n ma ttcr = 0.27 and f2yv = 0.73). The source sizes are plotted at the absorption redshift. The 
arrows represent continuum source upper limits, the open symbols are the 21-cm non-detections and the solid symbols are the detections. Right: The frequency 
of various source sizes, where the bold histogram represents the detections (upper limits are excluded). 



All 21-cm absorption detectio ns occu r in lower redshift 
(2 abs < 2.04) DLAs (see Fig.|4} and KC03 attribute this to these 
having a mix of both low and high spin temperatures (Fig.Q top), 
while the higher redshift DLAs generally have high spin tempera- 
tures. This may imply that the more distant DLA sample consists 
of a particular class of object. However, due to their large distances, 
none of the high redshift 21-cm sample have hosts which have 
been optically identified and indeed only one of the z < 2 non- 
detections has an optical identification (Fig.Q bottom). Therefore, 
we have insufficient information with which to determine the DLA 
host morphology and thus the likely covering factor. From this and 
the above points, we cannot rule out that perhaps the inferred (high) 
spin temperatures are an artifact of an overestimate of the covering 
factor. 

One way to address this issue is to estimate how the expected 
size of the intervening absorber compares to that of the back- 
ground continuum: In Table |2| we show the results from high res- 
olution radio images for the sources searched 5 and in Fig. [4] we 

5 Note that the continuum source sizes estimated by KC03 from the spec- 
tral energy distributions (SEDs) agree qualitatively with the results in the 



show the radio source size/absorption redshift distribution com- 
pared with various absorber sizes, D c . Although we are unc ertain 
which of these apply to a "typical" DLA (e.g. Brig as et alJll989t 
IChurchill & Charltonl[l999l: Ichengalur & Kanekarll2002l) . the bi- 
modal distribution of the continuum source sizes (particularly for 
the detections) itself may be indicative of a strong covering factor 
influence: Using the absorber host optical identifications for this 
sample we see a clear correlation between galaxy type and radio 
source size, where the spirals trace the upper segment (9qso > 1") 
of the distribution and the compact galaxies appear to be associated 
with the smaller radio sources (#qso < 0.2"). The fact that ab- 
sorption is detected only towards the extended continuum sources 
when occulted by a large galaxy strongly suggests that the cover- 
ing factor is a key issue. Toward the smaller continuum sources we 



table. Specifically, the flat SEDs for 0215+015, 0335-122 & 2342+342 sug- 
gest compact sources. For the remainder of their sources they use the ratio 
of the unresolved component's flux (Section \2. l\ to estimate the covering 
factor and, in order to avoid erroneous spin temperature estimates, exclude 
from further analysis the two sources (0957+561 & 1354+258) where this 
is expected to be low. 
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Table 2. Radio source sizes (arc-seconds) and VLBI and VLA morphologies of the QSOs illuminating DLAs searched for 21-cm absorption. 



QSO 


0QSO 


Ref 


VLBI 


Ref 


VLA 


Ref 


Comments [Ref] 


0235+164 


7.5 


16 


CJ 


25 


T 


16 


VLBI (5 GHz)/VLA (2 GHz) triple + extra strong component 


0248+430 


0.012 


22 


T 


22 


U 


1 


VLBI (1.6 GHz) weak core + 2 strong lobes or a CJ structure, 
outer J stronger than inner/VLA (5 GHz) < 0.8" 


0458-020 


3.6 


38 


CJ 


30 


CDom 


38 


VLBI (5 GHz)/VLA (5 GHz) core + weaker lobes 
(2 component, 1.4 GHz VLA [11]) 


0738+313 


0.009 


29 


CJ/T 


29 


T 


16 


VLBI (5 GHz)/VLA (1.6 GHz) triple, FRI like (core strongest) 


0809+483 


5.27 


3 






D 


3 


VLA (5 GHz) FRII like (strong lobes), as MERLIN 408 MHz [23] 


0827+243 


8.00 


17 


CJ 


39 


D 


17 


VLBA (2 & 8 GHz)/VLA (1.4 GHz) FRII like 


0952+179 


0.021 


39 


T 


39 






VLBA (2 & 8 GHz) either CJ or C-J-CJ (one side stronger) 


1127-145 


0.013 


24 


CJ 


24 


U 


1 


VLBI/VLA(1.4&5GHz) 


1157+014 


< 1.2 


6 






u 


6 


VLA (5 GHz) 


1229-021 


17.76 


4 






T 


4 


VLA standard triple, C + 2 weaker lobes 


1243-072 


10.17 


5 


CJ 


39 


T 


5 


VLBA (2 & 8 GHz)/VLA standard triple, core is dominant 


1328+307 


2.57 


14 


CJ 


39 


CDom 


14 


VLBA (2 & 8 GHz)/VLA (1.5 GHz) strong core + weaker component 


1331+170 


< 0.04 


9,31 


CJ 


41 


U 


9 


VLBA (2 & 8 GHz)/VLA (5GHz) 


1413+135* 


0.040 


20 


CJ 


20 


u 


1 


VLBI CJ + 2 more extended components on each side, core stronger 
VLA (5 GHz) 


1629+120 


2.025 


15 


D 


32 


T 


15 


VLBI (2 GHz) EVN/VLA 3 strong components, possibly FRII 


2351+456 


0.03 


22 


CJ 


22 


U 


1 


VLBI (1.7 GHz) strong 20-30 mas core + one weaker 100 mas jet 



0118-272 


2.0 


1 


T/mul 


34 


T/mul 


2 


0201+113 


< 0.007 


32 


C 


32 


CDom 


12 


0215+015 


9.45 


7 


C 


27 


T 


7 


0335-122 


< 0.4 


11 


C 


41 


U 


11 


0336-017 


8.0 


37 


CJ 


18 


T 


37 


0432-440* 














0438^136 


0.039 


35 


D 


35 


U 


1 


0439^133 














0454+039 


< 0.4 


11 


CJ 


41 


U 


11 


0528-250 


0.010 


40 


C 


40 


u 


1 


0537-286 


0.005 


40 


CJ 


40 






0906+430 


8.3 


14 


CJ 


22 


T 


14 


0957+561 


11.36 


28 


CJ 


21 


T? 


28 


1225+317 


< 0.24 


8 


C 


41 


U 


8 


1228-113* 


0.01 


41 


CJ 


41 






1354-107 














1354+258 


12.24 


15 






T 


15 


1451-375 


63.79 


33 


CJ 


39 


T 


33 


2128-123 


0.010 


26 


CJ 


26 


U 


1 


2223-052 


0.5 


13 


CJ 


40 


u 


4 


2342+342 


< 0.24 


36 






u 


36 



VLBI (5 GHz) core + extended 

VLA (5 GHz) at least 3 components, complex 

Core + weak extended (1% of total 21-cm flux within 2" [12]) 

VLBA (2 & 8 GHz)/ VLA (5 GHz) FRI, core dominated 

VLBI (2 GHz)/VLA (1.4 GHz) 

VLBI (5 GHz) unresolved/VLA (5 GHz) FRI, core dominated 

VLBI (5 GHz)/VLA (1 & 5 GHz) 

VLBA (2 & 8 GHz)/VLA (1 & 5 GHz) 
VLBI (5 GHz)/VLA (1 & 5 GHz) 
VLBA (2 & 8 GHz) 

VLBI (1.6 GHz) core strongest/VLA (1.4 GHz) core strongest 
VLBI (5 GHz)/VLA (5 GHz) lensed, looks like CLL core strongest 
VLA (2 GHz) possibly CJ/VLA (5 GHz) 
VLBA (2 & 8 GHz)/VLA: no map, 
may have flux on larger than VLBI scales 

VLA (5 GHz) FRI like, core dominates 

VLBA (5 GHz)/VLA (1.4 GHz) FRI like, core dominates 

VLBA (2 & 8 GHz)/VLA (1 & 5 GHz) 

VLBA (2 & 8 GHz)/VLA (1.4 GHz) 

VLA (8 GHz) 



Key VLBI morphology: CJ-core-jet morphology (most flux in core), D-double (2 lobes, no core visible), C-core, T-triple (core plus two 
lobes), mul = multiple components (> 3). N.B. it is not always possible to determine the core in the VLBI images. 

Key VLA morphology: U-unresolved, T-triple, usually the morphology is either FRI like (most flux in core) or FRII like (most flux in 
lobes), D-double, usually FRII like (core may be too faint), CDom-core dominated, usually a very strong core plus some fainter emission 
(some may look like very weak lobes). *Note that these are included in order to provide a comprehensive list but not included in the analysis 
(Fig.|4| - 1413+135 is an associated syst em and the lim its for 0432^140 and 1228-1 13 are poor (Appendix A) 

Re ferences: (lllUlvestad et alj<198lh. (2l|PerlevHl982l). ffllSchilizzi et al]ll982t) (4l lHintzen etai]|l983l) (5 fcower & HutchingsUl984l) . 
(6llStocke et all <1984l). m lAntonucci & Ulvestad <1985l) . (8 jlRogora et alj<1987l) . ( 9llBarthel et aljll988t). ( lOl lBriggs et alj<198? ), (11 ) 
iNeff & HutchingJll990l). (12^ IStanghellini et alj ll99cTri3)lFeies et alJ<1992l )"(14l |van Breugel et alj<1992) . (15llLonsdale et alj 19931) . 
f!6^ iMurohv et all ll993T) . (171 IPrice et alJ<1993l) . (IS yiGurvits et all j 19941) . C191 Lister et all Il994l) . 720> IPerlman et al] Il994 1, (21) 
ICampbell et all <1995l). (22) [ Polatidis et all <1995t) (23llReid et all <1995l)!(24^lBondi et alj <1996l) . (25^IChu et all Il996t) . (26llFev et ail 
| 1996|) . C27l |Fev & Chariot 1 19971) . C28l |Harvanek et all <1997l) . f 29^ IStanghellini et all <1997t). f30^IShen et all <1997l). f3nlBrowne et al] 
ll998t). (32)lDallacasa et allll998l). (33^ISaikia et alj|l998l). (34llShen et alj|l998t). (35llTmgav et alj|l998l) Oolwilkinson et alj<1998l) . 
f37^ lReid et alHl999l) . f38l lBarthel et alil2000l) . f39 Wev & CharroJfeOOOh . C40l lFomalont et alJfcOOCI) . (4D lBeaslev et 31112002^ 



Hi 21 -cm Absorption in DLAs 7 



would also expect spirals (together with compact galaxies ) since 
the covering factor would also be large. This is not observed. lKC03l 
do note that some of the identifications may be uncertain but, con- 
sulting the literature (e.g. lRao et all2003T) . we find similar identifi- 
cations, although the morphology of the z a b s = 0.09 DLA towards 
0738+313 is unclear l Cohen 2001) (this has therefore been flagged 
as unknown in Figs.Q|4]and[5j- The lower segment of the low red- 
shift bimodal distribution therefore consists of 3 dwarf/LSBs and 
3 unknowns 6 . That is, in this group there are as many compact as 
there are unidentified morphologies, though the latter could well 
turn out to be larger galaxies. Regarding the bimodal distribution, 
the binary probability of finding 6 out of 6 spiral galaxies towards 
sources of angular extent #qso > 0.5" and 3 (or more) out of 4 
non-spirals at #qso < 0.5" is just 0.5 per cent. Thus, it would seem 
there is some evidence to suggest a non-negligible covering factor 
effect in the current sample. 

Furthermore, as mentioned above, all but one of the non- 
detections at low redshift have yet to be identified optically, which 
excludes us from estimating a probable coverage, as done for the 
majority of the detections. However, in Fig.|4]we note that the low 
redshift DLAs not detected in 21 -cm tend to occult the largest back- 
ground sources. In fact 6 of these are located in the upper segment 
of the bimodal distribution, compared to only one in the lower seg- 
ment 7 . This strongly suggests that the low redshift non-detections 
may be mainly due to inadequate coverage. 

However, we emphasise one possible caveat to the above ar- 
gument. Several of the larger (Cqso > 1") background sources 
have complex morphologies, and in fact for 6 of the non-detections 
(0215+015, 0336-017, 0906+430, 0957+561, 1354+258 & 1451- 
375) these may be core dominated (Table|5J. The fact that there is 
no information on the core flux at the redshifted 21 -cm frequency 8 
nor on the extent of the absorbing gas means that we cannot un- 
ambiguously attribute the non-detections purely to low covering 
factors, although the smallest continuum source size is 8" at ^ 1.4 
GHz. 

One way to estimate the extent of the absorber required to 
fully cover the source is via the Compton limit of the emitting 
region, which places a lower limit on the size of the emitting re- 
gion accordi ng to the u pper limit of 10 12 K for the brightness tem- 
perature (e.g. Rantakvro et al. 2003). Wolfe, Gawiser & Prochaska 
(2003) apply this to 0201+113, where <9 QS o < 7 mas (Table 
|2j, to obtain a minimum extent of 4.9 mas (or 35 pc at z a b s for 
H = 75 km s" 1 Mpc" 1 , flatter = 0.27 and fi A = 0.73). Ap- 
plying this to the rest of the 21-cm sample gives 0.7-16 mas (or 1- 
130 pc at Zabs), all of which are small in comparison to #qso- The 
results of Wolfe, Gawiser & Prochaska (2003) are interesting, how- 
ever, since the extent of 0201+1 13 set by the Compton limit is very 
close t o that set by the radio obser vations at 1.6 GHz (< 2.5 x 5.0 
mas 2 . Hodges et alll9 84: Stanahellini et al. 1990). From studies of 
the bolometric background radiation in DLAs, Wolfe, Prochaska 
& Gawiser (2003) find that these absorbers consist mainly of the 
cold neutral medium (CNM, where T ~ 150 K and n ~ 10 
cm -3 ) rather than the warm neutral medium (WNM, T ~ 8000 
K and n ~ 0.2 cm' 3 ). This supports our argument for lower 

6 Excluding the unknown at 0.1" & 2 abs = 0.779 (2351+456). The 0.1" 
continuum size applies to the total extent, altho ugh the strong core flux 
originates within a 0.03" iPolatidis et alll995l) placing this closer to the 
other sources in this segment. 

7 If applying the 9qso = 0.5" cut-off this becomes 7 versus 3. 

8 These should, however, still provide a re asonable measure of the lower 
frequency source size te.g. lFeies et alll992t) . 



covering factors and hence lower spin temperatures in DLAs, al- 
though b ased on the high covering factor towards 0738+313 (Sec- 
tion Lane et alj <200Cl) : iKanekar et alj <200ll) suggest a two 
phase gas model dominated by the WNM (~ 75%) in both ab- 
sorbers. _nvis_js_£onsistent with the WNM+CNM mix in nearby 
dwarfs jYoungetall200lLl2003l) . where such low metallicity sys- 
tems are expected to reflect conditions in t he early Universe, with 
the values of [0 /H|~ 0.02 to 0.2 solar (e.g. lLipovetskv et alll999l: 
iLee et all2002j) spanning the range of typical DLA meta llicities up 
to ~ 3.4 IProchaska et all2003tlCurran et all2004t) . However, 
even with the warm gas contribution, the mean harmonic mean spin 
temperatures of ~ 800 K towards 0738+313 still lie in the "cool 
DLA", low redshift regime (Figs.Q&|5] below). 

In 0201+113 Wolfe, Gawiser & Prochaska (2003) find the 
CNM density to be n « 6 cm -3 and, by assuming that the to- 
tal Hi column density comprises several CNM (21-cm absorbing) 
clouds of column densities N ~ 2 x 10 20 cm -2 , they obtain an 
estimate of < 10 pc for each CNM cloud. While this argument has 
the draw-back that the total H I column density in the CNM phase 
must be estimated, it leads to a range of possibilities for 0201+1 13, 
the two extremes of which are: 

(i) The < 1 pc 2 1 -cm absorbing clouds are aligned with the line- 
of-sight towards the ~ 40 pc (at z a bs) background continuum, so 
that the coverage is low but the total column density is high. 

(ii) The clouds are distributed with little overlap, giving high 
coverage but a low effective column density. 

Having / ~ 1 and ATcnm <C -/Vhi could account for the 21- 
cm non-detections, while permitting a large fraction of CNM gas 
(Wolfe, Gawiser & Prochaska 2003). This emphasises how equa- 
tion Q can at best only provide a column density weighted mean 
harmonic spin temperature from a complex of clouds of various 
column densities, velocity dispersions and covering factors. 



3 CONCLUSIONS 

Upon a review of the literature, we suggest that less-than-complete 
covering factors may be at least partly responsible for the non- 
detection of 21-cm absorption in many DLAs. Although it is as im- 
portant as the spin temperature in determining the column density 
of the absorber, the covering factor is at best estimated from insuf- 
ficient data, and, more often than not, simply assumed to be unity. 
This assumption is prevalent in the high redshift cases, where high 
spin temperatures are claimed to be dominant, despite their being 
no DLA host identifications at z a b s ~ 0.9. 

Although all detections of 21-cm absorption occur at z a b s ?S 
2.04, nearly half of those searched at these redshifts remain unde- 
tected, thus giving an apparent mix of both low and high spin tem- 
peratures for this sample. However, of these non-detections, « 70 
per cent of the absorbers occult large background radio sources 
1"), thus increasing the likelihood of / < 1. Furthermore, 21- 
cm absorption tends to be detected: 

(i) Towards large background sources (> 1") only when the 
DLA is known to be associated with a large galaxy. 

(ii) Towards compact sources (< 0.1") when the DLA is asso- 
ciated with an LSB or dwarf galaxy. 

Naturally, we would also expect detections where the absorber is 
associated with a large spiral occulting a compact background con- 
tinuum, thereby providing large coverage. However, this is not ob- 
served. Of the 6 detections towards small background sources, 3 of 



8 S. J. Curran et al. 



5£, if) 



o 

O (N 



1 1 1 1 1 — 

- ▲ 


i i i i i 


-i 1 


1 1 1 1 1 




t 




t : 








t t ; 


:" t ft f 


t 




t * : 
t . 


★ 


t 






:t t 

★ 




t 


t 

t 


★ 

■ i ■ ■ I 





■ ' 





interpretation of low spin temperatures in these objects ( KCO^), but 
also with high coverage by spirals, a possibility strongly suggested 
by the work presented here. 

The lowering of covering factor estimates would introduce 
more uniformity to the spin temperatures of DLAs, while still sup- 
porting the hypothesis that large galaxies dominate the low red- 
shift population ( Ba ker et alj l2000l) with d warf galaxies consti- 
tuting more of the high red shift population lLanfranchi & Friacal 
2003 ). That is, unlike KC03 who suggest that this is evident through 
increased spin temperatures at high redshift, we hypothesize that 
the absorption cross-sections of compact galaxies at these redshifts 
simply lack the extent to effectively cover the background radio 
continua. 



Figure 5. Spin temperature/covering factor ratio against absorption redshift 
for the optically thin (r ^ 0.3) 21-cm searches. Again the FWHM of the 
non-detections is assumed to be 20 km s~ 1 (Section l2.2l . 



the hosts are unidentified, thus not ruling out the possibility of an 
approximately equal mix of large and compact DLA hosts. 

From the observed bolometric background, Wolfe, Gawiser & 
Prochaska (2003) propose that the 21-cm absorbing CNM is a ma- 
jor component of DLAs. They suggest that this is comprised of 
several lower column density clouds which may exhibit a large to- 
tal column density but low covering factor when aligned along the 
line-of-sight or, conversely, low individual cloud column densities 
and a large covering factor, where overlap along the line-of-sight 
is minimal. Therefore, when attributing the non-detection of 21-cm 
absorption in DLAs to the effects of the covering factor (or spin 
temperature), we should bear in mind that the total H I column den- 
sity may not provide a reliable estimate of the column density of the 
CNM along the line-of-sight to the background continuum. Further 
complications arise from the fact that insufficient detail is known 
about the flux contributions from the various components of the 
background continuum illuminating the DLA. The high resolution 
maps are at different frequencies than that of redshifted 21-cm and, 
if the components are extended, spectral indices may vary across 
them, making estimates of the effective continuum size at the ap- 
propriate frequency difficult. It is clear that high resolution obser- 
vations at the appropriate frequencies are required in order to fully 
address this issue. Combining these with the absorption profiles at 
variou s locations over the background continuum (as in lLane et alJ 
2000) will give far better estimates of the covering factor. 

Despite the above caveats, our analysis suggests that the cov- 
ering factor may play an important role, thereby significantly low- 
ering the spin temperature estimates of the DLAs for which the cov- 
ering factors have been assumed (i.e. most of the non-detections) 
and perhaps it would be more prudent to use T S p\ n f f versus absorp- 
tion redshift (Fig.[5} rather than just T sp i n , for an unknown value of 
/. The figure does in fact show that the values of T sp iu/ f span sim- 
ilar ranges for both the detections and non-detections with the spi- 
rals (stars) being more concentrated towards the low T sp iu/ f (and 
2aba) regime 9 (although, as stated previously, the majority of host 
identifications are unknown). This is apparently consistent with the 



9 Note however that the spread along the ordinate is considerably wider 
than in Fig.0(top): Over 1.2 dex in T sp i n // cf. 0.3 in T sp i n . Arguably this 
is mainl y influenced by 1229-021 at [0.39, 3.32] but also applies when we, 
as lKC03l include the associated system 1413+135 at [0.25, 3.28]. 
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Table 3. The known southern radio illuminated DLAs and sub-DLAs 
which fall into the Parkes 70-cm receiver band (B 1950.0 na mes are us ed 
throughout thi s pape r). T hese are obtained fro m iTurnshek et alj ll989l) a , 
lLanzetta et alj tl99ll) & and lEllison et alj l200ll) c . where A^hi and z a b a are 
the column density (cm -2 ) and the redshift of the DLA, respectively. v ^ B 
is the observed frequency of the redshifted 2 1 -cm line (MHz) and S is the 
flux d ensity of the background QSO at this frequency in Jy£see|Curraiy2^j3 
120021) . For 0528-2505 the value is the flux at 372 MHz ICarilli et all 19961) 
and for 1228-113 the value is from a 327 MHz Giant Metrewave Radio 
Telescope (GMRT) observation (Nissim Kanekar, private communication). 
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APPENDIX A 

Our initial motivation for this work was a search for 21 -cm absorp- 
tion in DLAs with the Parkes radio telescope. From the September 
2001 version of our catalogue of all known DLAs ICurran et alj 
2002) 10 we shortlisted those which are illuminated by radio-loud 
QSOs (i.e. those with a measured radio flux density > 0.1 Jy). From 
these, the DLAs with redshifts appropriate to the Parkes 70-cm re- 
ceiver bandwidth (i.e. 2.1 < z < 2.4) known to occult southern 
(S < 30°) QSOs were selected (Table 0. Although there are 7 
possible targets, we used the 90 hours scheduled to concentrate on 
3 of the 4 confirmed DLAs; 0432-440, 0438-436 & 1228-1 13 11 . 
0528-250 was excluded due to its relatively low flux density. The 
redshifts of these DLAs are known to within Az ~ ±0.002 (Sara 
Ellison, private communication), corresponding to Av ~ ±0.25 
MHz (w ±200 kms" 1 ). 

The observations were undertaken in January 2002 using the 
70-cm receiver on the Parkes 64-m antenna. We used the AT 
conversion system with the multibeam correlator, giving a band- 
width of 8 MHz (Az ps ±0.03) over 2048 channels (i.e. a spec- 
tral dispersion of 2.8 kms" 1 per channel) and position switch- 
ing with 1° beam throw. The weather was clear and dry giving 
system temperatures typically between 60 and 150 K. Although 



we restricted our observations to night time, we experienced se- 
vere radio-frequency interference (RFI) in this band. The RFI con- 
sisted primarily of sharp (FWHM < 6 kms) 'spikes' and broad 
(FWHM > 100 kms" 1 ) 'humps'. The humps in particular were 
found to shift in frequency between 120-s scans of the on- and 
off-source positions. Therefore, attempts to subtract the off-source 
scans from their corresponding on-source scans merely introduced 
strong RFI variations in the residual spectrum with characteristic 
length scales ~ 30 kms" 1 , i.e. the very scales typical of Hi 21-cm 
absorption lines in DLAs (see Table 2). The stronger spikes were 
also often accompanied by a powerful 'ringing' such that strong 
correlations between neighbouring channels could be seen to die 
away over ~ 300 kms" 1 scales. The data reduction was therefore 
carried out using software written specifically for the purpose of 
dealing with this RFI. The reduction steps are summarised as fol- 
lows: 

(i) The off-source scans were discarded and the on-source scans 
were low-pass filtered to provide a basic continuum shape which 
included the humps but ignored the spikes. The filtering scale was 
set so that ~ 50 km s" 1 absorption features could be recovered reli- 
ably. This continuum was then subtracted from the on-source scans 
leaving relatively flat residual spectra contaminated by the spikes. 

(ii) For each QSO, all > 500 contributing scans were corrected 
to the heliocentric frame and re-dispersed to a common frequency 
scale with a dispersion of 10 km s~ . 

(iii) A 1 a error array was constructed for each re-dispersed scan 
based on the rms in a sliding 20-channel window and channels with 
fluxes deviating by more than 3a were flagged and ignored in fur- 
ther reduction steps. 

(iv) Finally, the fluxes in corresponding channels in each scan 
were median filtered, values more than 2 a from the median re- 
jected and the remaining values combined to form the weighted 
mean flux density and 1 a error arrays. 

In Fig. [6] we show the resulting spectra for each DLA. It is 
clear that there are no significant, convincing absorption features 
over the mitigated intervals covered by the spectra (which include 
the frequencies expected from the optical redshifts). Generally, the 
noise level is quite high, w 0.1-0.2 Jy (cf. the theoretical value of 
a < 3 mjy per 10 km s" 1 channel) and some regions of the spec- 
tra were so severely affected by RFI that the errors are higher still. 
Where no reliable data could be gathered from any of the contribut- 
ing scans, we have assigned a normalized flux density of zero and 
an effectively infinite error. Note also that the low-pass filter ap- 
plied in stage (i) of the reduction procedure effectively prevents us 
from detecting absorption features with FWHM > 100 km s. The 
3 a optical depth limits quoted in Table Q for these sources must 
therefore carry these caveats. 

Due to the severe RFI, no flux densities could be obtained for 
our sources. Using the values in Table [3] for a — 0.2 Jy at 10 
kms" 1 resolution we obtain no 3a detection of the flux for ei- 
ther 0432-440 nor 1228-113 and r < 0.1 for the optical depth of 
21-cm absorption in the DLA towards 0438^-36. From equation|2| 
and the column density of the Lyman-a line (Table 0, this gives 
^£52. x FWHM > 3500 K km s~\ i.e. > 1200 K per 3 

km s" 1 channel or T s > 200 K for a FWHM of 20 km s" 1 (see 
Sectionl2~2l. 



10 A version of this catalogue is continually updated on-line and is avail- 
able from http://www.phys.unsw.edu.au/~sjc/dla 

11 1017+109, 1021-006 and 2136+141 are candidate DLAs which have 
yet to be confirmed using high resolution optical spectroscopy. 
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Figure 6. Parkes 21-cm spectra towards 0432^140 (top), 0438^136 (mid- 
dle) and 1228-1 13 (bottom). An uncertainty of Az ±0.002 corresponds 
to the ranges 430.6-431.1, 424.1^124.6 and 444.6^145.1 MHz, respec- 
tively, indicating that spikes have been successfully removed over the ranges 
expected from the optical absorption. 



